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Internal Josephson Effects in Spinor Dipolar Bose— Einstein Condensates 
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We theoretically study the internal Josephson effect, which is driven by spin exchange interactions 
and magnetic dipole-dipole interactions, in a three- level system for spin-1 Bose-Einstein condensates, 
obtaining novel spin dynamics. We introduce single spatial mode approximations into the Gross- 
Pitaevskii equations and derive the Josephson type equations, which are analogous to tunneling 
currents through three junctions between three superconductors. From an analogy with two inter- 
acting nonrigid pendulums, we identify unique varied oscillational modes, called the 0-tt, 0-running, 
running-running, 2mr & running-2'K, single nonrigid pendulum, and two rigid pendulums phase 
modes. These Josephson modes in the three states are expected to be found in real atomic Bose gas 
systems. 

PACS numbers: 03.75.Mn, 03.75.Kk 



I. INTRODUCTION 

The Josephson effect is a universal quantum phe- 
nomenon, defined as a current between two or more 
macroscopic quantum states with weak coupling driven 
by the relative phases of the macroscopic wave functions. 
The effect for a junction of two superconductors was orig- 
inally predicted by Josephson [l[ and discovered by An- 
derson and Rowell Q. When the effect was first pre- 
dicted and discovered it was considered a characteristic 
phenomenon of superconductors. However, the effect has 
been found in diverse fields since Feynman redefined the 
Josephson effect as the quantum tunneling between two 
levels 0|. In particular, Maki and Tsuneto have consid- 
ered nonlinear ringing between equal spin pairing states 
| ||) and | ||) in superfluid 3 He-A as a Josephson junc- 
tion between the internal degrees of freedom [J]. Webb 
et al. subsequently observed the effect by NMR Q . 

Such Josephson effects have also been studied in atomic 
Bose-Einstein condensates (BECs). Smerzi et al. derived 
Josephson type equations in a double well potential from 
Gross-Pitaevskii (GP) equations with two-mode approx- 
imations, revealing three oscillations, namely the 0, tt 
and running phase modes |f| 0] • Macroscopic quantum 
self trapping from the running mode and AC Joseph- 
son effects from the phase mode have been observed 
by Albiez et al. [8[. Furthermore, Zhang et al. consid- 
ered transitions between three energy states for spin-1 
BECs driven by spin exchange interactions as internal 
Josephson effects |9j and Chang et al. observed the inter- 
nal Josephson oscillation ITo| . 

In addition, studies have focused on the magnetic 
dipole-dipole interaction (MDDI) in BECs. The inter- 
action between spins has a characteristic symmetry for 
spin and orbit. Therefore, the interaction is theoretically 
expected to give new quantum phases [ll|, [l2|, [l3j |. the 
Einstein-de Haas effect [14[ and so on. Griesmaier et 



al. realized spinor dipolar condensates using 52 Cr atoms 
[l5[ which have larger magnetic moments than alkali 
atoms. The condensates clearly show shape-dependent 
anisotropy of the MDDI [H, E3 • Thus, the study of the 
MDDI has opened new avenues for investigating spinor 
condensates. 

For the MDDI, we can demonstrate new internal 
Josephson effects in spin-1 BECs. The internal Joseph- 
son effect without the MDDI is analogous with two junc- 
tions between three superconductors because the weak 
coupling, namely the spin exchange interaction, permits 
transitions only from |0) to | ± 1) and from | ±1) to |0). 
However, it is possible to transfer from | ± 1) to | =p 1) if 
the MDDI drives the transitions because the anisotropic 
interaction breaks the conservation law of magnetization, 
as will be discussed later. As the previous study, Cheng 
et. al. have discussed the Josephson effects driven by 
the MDDI in the conservation law [18[- Therefore, we 
consider out of the law that the transitions through the 
MDDI are analogous with Josephson effects for three cir- 
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FIG. 1: (color online) The transitions, driven by the Hamilto- 
nians of a spin exchange interaction Hse (black arrows) and 
an MDDI Had (red arrows), between the energy states (solid 
line) are analogous to the Josephson effect for cooper pairs 
tunneling through (a) two and (b) three insulators (white 
boxes) between super conductors, named SCI, SCO, and SC-1 
(red boxes), (a) and (b) compare differently with the MDDI, 
was the MDDI can also drives transition from | dz 1) to | =p 1). 
The blue circles represent the condensate. 
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cular junctions between three superconductors (see Fig. the equation of the spinor |£) = (£1, £-i) T in H = Hz, 

®- . 

(3) 



In this paper, we introduce the spin dynamics in spinor 
dipolar BECs in Sec. [Ill using the GP equations with 
single mode approximations, and we demonstrate new 
Josephson oscillations in Sec. IIII1 Section HVl is devoted 
to a conclusion. 



II. FORMULATION 

A. Single mode approximation 

Before we consider the internal Josephson effect for 
spin-1 BECs, we analyze the spin-1 GP equations with 
an MDDI: 

a + C2-T • F a 

where the integers a and (3 represent the spin Oand ± 1 
states, V(r) is the trapping potential, and \i the chemical 
potential. The spin density vector F = (F Xl F y , F z ) and 
the vector of the spin components F a p = (F£p, F^p, F^p) 
are represented by F i=Xy y^ z = i/j^F^pipp given by the com- 
ponents F^p of the spin matrices Fi for spin-1. The Zee- 
man term is presented by Lande's ^-factor of an electron 
g e: Bohr magneton /ig, and external magnetic field H. 
The coefficients Co = (go + 2^)/3 and = (g2 — 9o)/3 
are the interaction parameters for gi = 47rh 2 ai/M rep- 
resented by the s-wave scattering lengths a\. The coef- 
ficient of the MDDI is given by Cdd = PoglPs/^ w ^ n 
e = (r — r')/|r — r'| being a unit vector. 

Since it is difficult to derive the internal Josephson 
effect from Eq. ([1]) directly, we introduce the single mode 
approximation [9[: 



^(r,t) = V7V^(t)0(r)exp 



ipt 



(2) 



where <fi satisfies the eigenvalue equation (—h 2 V 2 /2M + 
V + conjcf) = fKp with the relation j dr\(j)\ 2 = 1. The ap- 
proximation can be used when the shapes of the conden- 
sates are decided by the spin independent terms, namely 
M > |c2 1- The condition is satisfied for 87 Rb and 23 Na. 
The validity of the approximation in spin-1 dipolar BECs 
has been discussed by Yi and Pu [19]. They found that 
the approximation is affected by the anisotropic param- 
eter of the shape of the condensates A and the ratio 
Cdd/ 1 02 1- For small A, the approximation is well defined 
even if the ratio is large. Here, we assume that the ratio 
is smaller than 1, allowing the approximation to be used. 
Hence, introducing Eq. ([2]) into Eq. (pQ), we can derive 



ifrj t \0 = (Hz + H SE + H dd )\0, 



given by the Zeeman, the spin exchange interaction, and 
the MDDI Hamiltionians, 

( H \ 

Hz = -gfJLB (4) 
\ -H ) 

1 - 2p_i p_i 
Hse = c ( po-i 1 - Po Poi I (5) 
pio 1- 2pi 

Po P-io 
Hdd = Cdd I Po-i 1 - Po Poi 
pio Po 

pio Po 
+ c dd I Poi P-n+Pi-i Po-i 
Po P-io 

/ rn \ 
+2c z dd 0, (6) 
\0 -m J 

where pij = (pa = pi) are components of the density 
matrix p = |, c± = c x dd ± c y dd with c z dd given by 



i °dd AT 

c d d = -yJV 



r'|3 



(7) 

c = C2N J dr\(j)\ 4 are the interaction parameters, and m = 
pi — p-i is the magnetization. Comparing Eq. ([5]) with 
Eq. (|6]), we can conclude that the term with c dd in H dd 
includes the operators |zbl)(=pl|, which project from |=pl) 
to |±1), where |1) = (1,0, 0) T , | - 1) = (0,0, 1) T are unit 
vectors. These operators are not included in Hse- We 
again emphasize that the transition from | ± 1) to | =p 1) 
is possible only by the MDDI. 

Using the relation ^ pi = 1 we can rewrite the spin 
exchange and the dipole term to give 



Hsi 



dd 



Po P-10 
( c + c dd) I Po-i 1 - Po Poi 
pio po 




(8) 



Equation ((5)) clearly shows that a c term can be included 
in the c dd and c dd terms. In short, the MDDI is partly 
renormalized in the spin exchange interaction. 

In conclusion for the MDDI, three important roles for 
spin dynamics appear in Eq. (|6|). The first is the spin 
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exchange effect, as in Eq. (j5j), given by the c dd and c dd 
terms of Hdd- The second is the transition from the ±1 
to the =pl states, given by the c dd term. The last is an in- 
teraction between the magnetization and the spin, given 
by the c dd term, which implies that the magnetization 
produces a molecular field. 
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B. Josephson type equations 

In order to investigate the Josephson effect for tran- 
sitions between the three states, we substitute £j = 
y J~Pje %9j in Eq. ([3j), deriving the Josephson type equa- 
tions for po, m and relative phases 9 = 9\ + 0-\ — 29 o 
and rn — 9\ — 



FIG. 2: (color online) Analogy between the dynamics of 
Josephson junctions an d two interactin g nonrigid pendulums 
with lengths Iq oc d po\/ (1 — Po) 2 — m 2 and Z# m oc (c^ d /c)le. 



whose Hamiltonian is represented as 



2 

Po -c'p y/(l - po) 2 -m 2 sin9 

2 

+ ^ c ddA){(l - Po) sin<9cos<9 m + racos<9sin# m }, (9) 



H = Hq + He n 



H 



int i 



(13) 



m -^ddPoWi 1 - Po) 2 ~ m<2 sin6> m 
+ (1 — po) cos#sin# m + msin#cos# m }, 



(10) 



2 , n 9 v , 2 / (l-po)(l-2p )-m 2 

= -c (1 — 2p ) + t c / — cost^ 

h h - po) 2 - m 2 



h° dd 



(l-p )(l-2p )-m 2 



COSt 



po) 2 ~ m 2 
(1 — 2po) cos 9 cos 9 m — m sin 9 sin 9 m > , 



2/^ . „ v . 2 , rap 
--(Ez + c m) + -c — . 



(11) 



— c 



mpo 



dd 



Vi 1 - Po) 2 -m' 



: COSt 



po sin # sin t 



(12) 



where ^ = —g e pBH is the Zeeman energy and c' = 
c + c dd and c" = c + 2c^ d are the interaction coefficients. 
The equations clearly demonstrate the definition of the 
Josephson effect, namely a current between states driven 
by a relative phase. 

Firstly, from Eq. ([TOj) we can easily see that conser- 
vation of the magnetization is broken by the c dd term. 
When c dd = 0, the magnetization becomes a conserved 
value. Under the single mode approximation, c dd — for 
uniform density and a spherical shape, which is obtained 
from the integration, Eq. ([7j). Secondly, we note that the 
variables po and m are canonical conjugates of and m , 
respectively. Therefore, these equations can be derived 
from the canonical equations of motions: 

• - _ 2 dH e- 2dH m - 2 dH and 9 - _ 2 dH 

^° h 80 ' h dpo ' m h 00m ' m h dm ' 



where 
He 

Ho rn : 



c'{p (l - Po) + PoVi 1 - Po) 2 - m 2 cos 9}, (14) 



—m 2 + £ 2 ra + c dd po^f (1 - po) 2 - ?™ 2 cos0 m , 

(15) 

Hint = c^{p (l - Po) cos6>cos6> m - p ^sin6>sin(9 m }. 

(16) 

Since the Hamiltonian H(po(t),m(t),9(t),9 m (t)) is not 
explicitly time dependent, the relation dH/dt = 
dTL/dt = is satisfied. A comparison of the Hamil- 
tonian with that of a pendulum of length I and an- 
gle ip, H p (P,(p) = P 2 /2M + Mgl cos ip, shows that the 
dynamics described by the equations may be equal to 
that of two nonrigid pendulums, because He and He m 
are similar to the Hamiltonian of nonrigid pendulums 
with lengths proportional to c! po\J (1 — Po) 2 — and 
c dd po\/ (1 — Po) 2 — 7T^ 2 and angles 9 and m , respectively. 
The two pendulums have some interaction with each 
other, represented by terms in Hint that include the prod- 
uct of 9 by 9 m . Therefore, we consider that Eq. ([T3]) 
represents the Hamiltonian of two interacting nonrigid 
pendulums (see Fig. [2|). 



III. SPIN DYNAMICS 

In this section, we discuss solutions of the equations 
for Ez — 0, namely small oscillations of two nonrigid 
pendulums in Sec. IIII A\ oscillations of a single nonrigid 
pendulum in Sec. IIII Bl and oscillations of two interacting 
nonrigid pendulums in Sec. IIII CI 



A. Small oscillation around stationary solutions 

Not being able to solve the Josephson type equation 
directly, we first obtain the solutions for the stationary 



4 



states: 



po Pa, 



and 



Po 



1 

2' 



0, 6>=(2n + l)7r, 



rn = 0, = 2n7r, # ri 



n7r. 



(17) 



(18) 



Here p a and # a are arbitrary constants. For Ez = 0, 
we cannot obtain the stationary magnetization m ^ 0. 
In order to study the features of the oscillation of the 
pendulums, we perform linear analysis by introducing the 
approximations po ~ pos+^Po? m — m s + Sm, c± S + S0 
and # m ~ # ms + where these deviations from the 
stationary solutions are very small, into the Josephson 
type equations. Linearizing the equations by ignoring 
squares of the deviations, we obtain the equations for 
the deviations; the solutions around Eq. (jTTJ) are 



Sm = SO = 
and those around Eq. (|I8j) are 



Spo 
Sm 



u Po± 



Spo 
Sm, 



(19) 



(20) 
(21) 



where u%± = 4(c' ± c dd ) 2 /h 2 and u 2 m± = 4c dd {Tc z dd ± 
(c dd ± c dd )/2}/h 2 with ± indicating the equations for 
# ms = 2n7r and (2n + l)7r, respectively. From Eq. (fT9]h 
we can conclude that there are no stable oscillational 
solutions around Eq. ([T7|) . On the other hand, Eqs. 
([20]) and ([2T]) are easily solved and have the solutions 
Spo = A p cosa;p±t and Sm = A m coso; m ±£, whose ampli- 
tudes are A p and A m . Considering the deviations from 
Eqs. ([T9]) to ([2T]h we can conclude that there are no oscil- 
lations around = (2n+l)7r, namely the tt phase mode in 
the canonically conjugate variable {po,9}, whereas there 
are and tt phase modes in {m, m }. 



B. Josephson effect for a nonrigid pendulum 

Next, in order to review the 0, 7r, and running phase 
modes in spin-1 BECs, we study the Josephson effect for 
c dd = obtained for spherically shaped condensates. Un- 
der this condition, the spin dynamics are represented by 
Eqs. ([9|) and ([TO]) with c dd = 0, where there is no need to 
consider m and m because m is constant. The equations 
are the same as for the Josephson effect driven by the spin 
exchange interaction in Ref. [9[ if the quadratic Zeeman 
effect, which is effective for the hyperfine interaction be- 
tween the nuclear and electron spins in the Zeeman splits 
fToL [20| , is introduced in the equation. Therefore, con- 
sidering the quadratic Zeeman effect Sdd cx H 2 ^ given 
by the effective dipole magnetic field H e ff cx c dd m, we 
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FIG. 3: (Color online) Contour lines of Eq. ([24]) for (a) m = 
and (b) 0.3 with Sdd = 0.3c' (> 0), represented as white. The 
blue, purple and magenta lines show the 0, running and ty 
phase modes, respectively. 



obtain 



Po 



- c' po V(I - po) 2 - rn 2 sin 

2 2 
= —^6 dd + -J(l-2p ) 



(22) 



2 (l-p d )(l-2 Pd )-m« 
h ^(1 - Po) 2 - m 2 

where po is the canonical conjugate of 0. Hence, these 
equations are derived from the canonical equations of mo- 
tions for the Hamiltonian: 



W = S dd (l- po)+c{po(l- po) 



+Poa/(i" 



Po 



) 2 — m 2 cosO. 



(24) 



The Hamiltonian directly represents the system of a non- 
rigid pendulum. Using W = H' = W(po(0), 0(0)) for 



TL' = 0, we can solve Eq. (|22|) and obtain the solu- 
tions expressed by Jacobi's elliptic function cn(a, fc); for 



Po = Pb 



+ (pc ~ Pb)cn 2 (^\f2c'8 dd (p c - ~pa)t, k^j , (25) 



and for c dd < 



Po = Pb - 



(pb-pa)cn 2 (^\/-2c'5 dd (p c - p a )t, k^j , (26) 



where p a ,b,c are the roots of po = 0, and k = {(p c — 
Pb)/(pc — Pa)} 1 ^ 2 - Pendulum- like oscillations occur on 
the contour lines of Eq. ([24]) for energy being conserved, 
shown in Fig. 03 The dynamics can be classified into 
three modes; the dynamics on the lines around = 
and 7r are called the and tt phase modes respectively 
and the line from = — 2tt to 2tt is the running phase 
mode. The phase mode corresponds to the motion of a 
pendulum oscillating around = with a varying length. 
The 7r phase mode also shows an oscillation around 6 = 
7r,a characteristic dynamics for the nonrigid pendulum. 
Finally, the running phase mode represents the rotational 
dynamics. 



C. Josephson effect for two pendulums 
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FIG. 4: (color online) Typical 0-tt phase mode. The dynam- 
ics of the four variables po, m, and m in (a) and (b) and the 
phase projections of the motions in {po, 0} and {ra, m } in (c) 
and (d) are given for p (0) = 0.51, m(0) = 0.01, 0(0) = tt/18, 
m (O) = 19tt/18. 



Here we discuss the dynamics in Eqs. (|9j) to ([T2]) . 
Solving the equations numerically using the fourth-order 
Runge-Kutta method, we obtain several characteristic 
results for two interacting nonrigid pendulums. In the 
calculations, we estimate that the interaction parame- 
ters satisfy the relations c f /c dd = 11 and c" /c dd = 12, 
which are determined from the order estimations c 



dd 



c dd ^ c dd an d c °° ^ C dd- ^ ne orcl ers between the spin ex- 
change interactions and MDDIs have been discussed by 
Yi and Pu [2l|. In this subsection, we denote the pen- 
dulums having angles 6 and 6 m as 6 and m pendulums, 
respectively, considering the dynamics from the time de- 
velopments of the canonically conjugate variables (in Sec. 
IIIIC ip and the Hamiltonian (in Sec. IIII C 2j) . 



1. Modes in the dynamics of four variables 

First, we present the typical 0-tt phase mode for which 
and 6 m oscillate around = and m = 7r, respectively; 
shown in Fig. 0] (a) and (b). The phase mode is clearly 
obtained by applying small deviations from the station- 
ary solution, Eq. ([T8]) . Considering the motion of the 
pendulums, the and m pendulums make small oscil- 
lations around and tt as the lengths grow and shrink. 
The features of the oscillations are represented by trajec- 
tories projected in the phase spaces {po, #} Fig. H](c) and 
{m,0 m } (d). Since the trajectory is confined to a simple 
closed orbit, we can see that the motion of the pendu- 
lums is similar to that of two independent pendulums, 
namely a periodic motion. On the other hand, changing 
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FIG. 5: (color online) Quasi-0-7r phase mode for po(0) 
0.51, m(0) = 0.36, 0(0) = tt/18, m (O) = 19tt/18. 



m(0) = 0.01 to m(0) = 0.36, we obtain the quasi-0-7r 
phase mode, shown in Fig. [5l The trajectory plotted in 
Fig. [5] (c) is not an orbit and hence we regard the spread 
trajectory as evidence of chaotic motion of the pendu- 
lum. In order to investigate the chaotic motions, we plot 
Poincare mappings in Fig. [6j The mappings represent 




FIG. 6: Poincare mappings for the 0-tt modes within 
(cj d /h)t = 500. The upper and lower figures show the projec- 
tions of the Poincare sections at T = 2ty/co po - for the phase 
spaces {po,0} and 2it/uJm- for {m, m } in Fig. |4]and Fig. El 
respectively. 



projections of the cross sections at t = nT (n = 1, 2, . . .), 
where the period T is given by eigenfrequencies in the sys- 
tems, providing direct evidence of chaotic motion for the 
0-7T phase modes, because there is only one point in the 
Poincare mappings where the motions are not chaotic but 
periodic. However, it is unexpected that the typical 0-tt 
phase mode should be chaotic. Therefore, we conclude 
that the motions of the two pendulums are essentially 
chaotic and thus the dynamics do not occur on a contour 
line of the phase mode. The m pendulum is also not 
a dynamic on a contour line of the tt phase mode. 

Setting m(0) = 0.37, where the other initial conditions 
are same as the 0-tt phase mode, we can produce the 0- 
running phase mode in Fig. The phase mode shows 
that the pendulum oscillates around 6 = with a small 
amplitude chaotically, whereas the m pendulum rotates 
continuously, given by the dynamics of m in Fig. [71 (b). 
The m pendulum has a long stay around m = 2n7r, 
which can be understood by Fig. [71 (b). 

Thus, we can obtain the transition from the 0-tt phase 
mode to the 0-running phase mode by changing m(0). 

Third, the running-running phase mode is shown in 
Fig. [HI The solutions show that the motions of the pen- 
dulums are rotations with different angular frequencies. 

Fourth, we present the 2nn & running-2i\ phase mode 
in Fig. where n is an integer. In this mode, the dy- 
namics of the two pendulums show strange and inter- 
esting motions. Especially, the pendulum repeats an 
oscillation around = 2nn and a rotation from = 2nn 
to (2n + 1)tt. The dynamics indicates a transition from 



FIG. 7: (color online) 0-running phase mode for po(0) = 
0.51, m(0) = 0.37, 0(0) = tt/18, m (O) = 19tt/18. 
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FIG. 8: (color online) running-running phase mode for 
po(0) = 0.6, m(0) = 0.2, 0(0) = tt/2 and m (O) = tt. 



the contour lines of the 2tt phase mode to another line of 
the running phase mode, which is shown clearly in Fig. 
[9](c). Due to this transition, transitions should occur for 
H = in {m,0 m }. Also, the m pendulum has a time 
average of (0) = 2ir. However, the pendulum repeats a 
rotation from ~ to 4tt and returns to ~ 0. 

We obtain the motions of a single nonrigid pendulum 
in Fig. [10] and two rigid pendulum in Fig. [TTJ namely 
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FIG. 9: (color online) 2n7r & running-2ix phase mode for 
po(0) = 0.1, m(0) = 0, 0(0) = tt/2, m (O) - tt/3. 
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FIG. 10: (color online) single nonrigid pendulum phase 
mode for p (0) = 0.2, m(0) = 0, 0(0) - tt/2, m (O) = tt. 



the single nonrigid pendulum and too rigid pendulums 
phase modes. For the single pendulum, the pendulum 
oscillates around = 0, whereas the length of the 9 m 
pendulum grows and shrinks without an oscillation be- 
cause m = 0. The phase space in Fig. [10] (c) shows a 
single trajectory, indicating dynamics like that of a rigid 
pendulum. On the other hand, Fig. [11] (a) and (b) show 
the running-running mode for rigid pendulums given by 
the constants po and m. Then, Fig. [11] (c) and (d) are 
very difficult solutions for understanding the dynamics. 




0.01 



FIG. 11: (color online) The motion of two rigid pendulums 
for the initial conditions po(0) = 0, m(0) = 0.8, 0(0) =0 and 
m (O) for (a) and (b), and p = 0, m = 0.01, tt/4, 
and 0m = tt for (c) and (d). 



The pendulum exponentially deviates from the initial 
angle to = tt. The m pendulum, however, exhibits the 
rotation of a rigid pendulum. 



2. Contour lines of the Hamiltonian 

Here we discuss the various modes for two pendulums 
from the Hamiltonians. 

It is impossible to plot the contour lines for the to- 
tal Hamiltonian, namely Eq. ([T3]h because of the four 
dimensions given by po? m i Q an d m . By considering 
the Hamiltonian of the pendulum Tie and the m pen- 
dulum TLe m1 however, the features of the motions of the 
pendulums can be obtained. 

As mentioned in Sec. Ill Bl TL is a conserved value. 
However, TLq^ He m and Hint are not conserved. Figure 
[T2l shows the time development of these Hamiltonians, 
indicating the non-conservation, noting though that the 
values in (a) exhibit only small oscillations that do sug- 
gest energy conservation. Therefore, the trajectory in the 
phase spaces does not have to follow a contour line. As 
examples, we show the trajectories and contour lines for 
the 0-7T and 2nn & running-2ix phase modes in Fig. [13] 
Naturally, the contour lines change with time because TLq 
and Ho m include the parameters m(t) and po(t), which 
are time dependent. In Fig. [131 we plot the lines for the 
time averages m = (m) and po = (po) in He and He m - 
The trajectories of the 0-tt phase mode almost follow the 
contour lines in Fig. [13] (a) and (b). On the other hand, 
for the 2nn & running-2ix phase mode of (c) and (d), the 
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FIG. 12: (color online) Time development of the Hamiltoni- 
ans. (a), (b), and (c) have the same initial conditions as for 
Fig. |U Fig. [5] and Fig. [9] respectively. 
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tion from the to running phase modes is clearly seen in 
{po,0}. An important result is that the dynamics of the 
m pendulum cannot be discussed as transitions between 
the 0, 7r and running phase modes. 



IV. CONCLUSION 

By introducing the single mode approximation into the 
spin-1 GP equations with the magnetic dipole-dipole in- 
teractions, we show the transitions between the three 
states, which are analogous with Josephson junctions of 
a ring of three superconductors. Furthermore, deriving 
the canonical equations of motion from the Josephson 
equations, the Josephson junctions are found to also be 
analogous with two nonrigid pendulums having interac- 
tions with each other. First, in order to consider a sim- 
ple nonrigid pendulum, we assume that the shapes of 
the condensates are spheres. In this way, we review well 
known three phase nodes, the 0, ir and running phase 
modes. Second, we numerically solve the equations for 
two nonrigid pendulums, showing several motions: the 
0-7T, 0-running , running-running , 2nir & running-2ir, 
single pendulum, and two rigid pendulum phase modes. 
Finally, we discuss the transition between the modes from 
the non-conserved Hamiltonians 7Y<9, 'He rri , and Hint- 

We consider that the Josephson effect, which has 
been found in various phenomena from condensed matter 
physics to classical physics, is a universal and important 
physical phenomenon. In this study, the effect was dis- 
cussed only in relation to BECs. However, we expect 
that the study of the Josephson effect of the three states 
will be useful in many other fields beyond atomic BECs. 



FIG. 13: (color online) Trajectories and contour lines for the 
0-7T, (a) and (b), and 2n7T & running-27r phase modes, (c) 
and (d). The trajectories are shown as green dots in the 0-tt 
mode and red dots in the 2n7r & running-2-K mode. The lines 
are given by Eq. ([T4]) for m = (m) ~ in (a) and (c) and 
Eq. dHJ) for p = (p ) ~ 0.5 in (b) and 0.47 in (d). The 
magnitude of the Hamiltonians increases from black to white. 



trajectories do not follow the lines. However, the transi- 
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